function run()
  global R1 R2;
  global t1 t2;
  
  tic
  R1 = h2(64, 0.16, 64, 1, 10000, 0.1, 0.1);
  t1 = toc;

  cd ..
  tic
  R2 = h2amz(64, 0.16, 64, 0.0001, 1, 0.1, 0.1);
  t2 = toc;


end
